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Abstract. In order to identify possible experimental signatures of the superfluid to Mott- 
insulator quantum phase transition we calculate the charge structure factor S{k,Lj) for the 
one-dimensional Bose-Hubbard model using the dynamical density-matrix renormalisation 
group (DDMRG) technique. Particularly we analyse the behaviour of S{k,(jj) by varying — 
at zero temperature — the Coulomb interaction strength within the first Mott lobe. For strong 
interactions, in the Mott-insulator phase, we demonstrate that the DDMRG results are well 
reproduced by a strong-coupling expansion, just as the quasi-particle dispersion. In the super- 
fiuid phase we determine the linear excitation spectrum near k = 0. In one dimension, the 
amplitude mode is absent which mean-field theory suggests for higher dimensions. 



Recent experimental realisations of optical lattices make it possible to investigate the 
properties of ultracold dilute atoms in a new regime of strong correlations (TJ [2]. Tuning 
the strength of the laser field, the effective interactions between the atoms can be tuned to 
be stronger than their kinetic energy. The competition between the Coulomb energy and the 
kinetic energy may even drive a quantum phase transition between superfluid (SF) and Mott 
insulating (MI) phases. The Bose-Hubbard model (BHM) captures the essential physics of this 
problem. In previous work [3] we studied the photoemission spectra of the one-dimensional (ID) 
BHM both in the MI and SF phases. In this contribution we extend our preceding investigations 
of the dynamical properties of the BHM by analysing the dynamical structure factor S{k,uj). 
Analytically, S{k, oj) has been determined by mean-field approaches [H [S] and, numerically, for 
rather small ID systems, by exact diagonalisation [U] or, at finite temperatures, by QMC 0. 
Here, we carry out large-scale dynamical density-matrix renormalisation group (DDMRG) [H] 
calculations in order to examine the density-density correlations of N interacting Bose particles 
on L sites at zero temperature. Moreover, fixing p = N/L = 1, we set up a perturbation theory 
whose results can be compared with the DDMRG data in the strong-coupling limiting case. 

The Hamiltonian of the ID BHM reads 

where b] and b- are the creation and annihilation operators for bosons on site j, hj = b\b, is 
the boson number operator on site j. The tunnel amplitude between nearest neighbour lattice 
sites is denoted by t, and U gives the on-site Coulomb repulsion of the Bose particles. For the 
DDMRG treatment of ([1]) we use periodic boundary conditions, take into account = 5 bosons 



per site and keep up to m = 500 density-matrix eigenstates, ensuring that the discarded weight 
is always smaller than 3 x 10"^. 

The dynamical charge structure factor is defined as 

S{k,Uj) = '^\{'ll^n\nk\ll^o)\'^6{u; - LOn) , (2) 
n 

where \ipo) and \ipn) denote the ground and n-th excited state, respectively, and ojn = En — Eq 
gives the corresponding excitation energy. S{k,uj) characterises the density-density response of 
the Bose gas and is directly accessible in experiment, for instance, by Bragg spectroscopy |9]. 
Prom lowest-order strong-coupling perturbation theory we find for p = N/L = 1 (N, L — t- oo) 



5(^) {k, a; > 0) = ' y " ^ sin2 kS{u;-U + 2t cos W5 + 4cosfc) . 

Then, for U - 2V5 + 4cosA; < u < U + 2V5 + 4cos/c, we have in units of t 
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where < A; < 2it. Since DDMRG provides S{k, oj) with a finite broadening r], it is 
useful to convolve the strong-coupling result S^^\k,ijj) with the Lorentz distribution |10| : 
Snik,uj) = J^^du'S{k,uj')r]/ [tt [{oj - io'f + 7]'^]]. For the ID BHM, from Eq. (g]), we thus 
get 



U + 2AV5 + 4cosA;)2 + r?2 ' 



(5) 



Figure [U illustrates the frequency and momentum dependencies of S{k,u}). Peaks in S{k,u) 
assign charge excitations. For A; ~ 0, S{k,uj) exhibits a maximum around u k, U that can be 
attributed to excitations across the Mott gap. Within strong-coupling theory @ this signature 
stays at uj K.U for all k but becomes more pronounced as k reaches the Brillouin zone boundary. 




Figure 1. Frequency dependence of the dynamical charge structure factor S{k,Lj) at various 
momenta k = 2Trj/L (for j = l,2,...,L/2 from bottom to top). DDMRG data (symbols) 
were obtained for L = 16; S^^\k,io) (lines) give the corresponding first-order strong-coupling 
expansion results (for the same value of r] = 4t). The dashed line indicates the point oo = U . 




Figure 2. 77-dependence of the dynamical structure factor S{k, lo). DDMRG data (symbols) are 
compared with the corresponding strong-coupling results (dashed lines) at various momenta k. 
The solid green line shows the unconvolved result of the strong-coupling expansion. 



Our DDMRG data corroborate this prediction (cf. the left panel of Fig. [T] depicting the results 
for t/U = 0.02). Naturally, as U gets smaller, differences between the numerical DDMRG and 
the analytical strong-coupling results emerge. Most importantly, the position of the maximum 
in S{k,uj) varies with k: it is shifted to smaller (larger) frequencies for k values near the band 
centre (Brillouin zone boundary); see middle and right-hand panels of Fig. [TJ We have to include 
higher-order t/^7-corrections to reproduce this feature analytically. Note that the maximum in 
S{k,uj) amplifies as A; — )■ vr. 

Next, we investigate the dependences of S{k, uj) on the broadening rj and the system size L to 
scrutinise whether the DDMRG data "converge" to the unconvolved strong-coupling result ^ as 
7] ^ and L — )• 00. Figure [2] — showing S{k,uj) in the MI phase at A; = vr/S (left), tt/2 (middle) 
and vr (right panel) for L = 16, 32, 64 r/ = 4i, 2t, t, respectively — demonstrates that this is 
indeed the case. Firstly, for t/U = 0.02, we see that the DDMRG results are in a satisfactory 
accord with si^\k,uj), where as a matter of course for smaller system sizes a larger value of rj 
has to be used to achieve good agreement. Secondly, S{k,uj) calculated by DDMRG approaches 
the curve given by Q for increasing system size L and decreasing broadening 77. 

Finally, we look at the changes in the dynamical density-density response as the system crosses 
the MI-SF quantum phase transition with decreasing Coulomb interaction strength. Figure [3] 
shows the intensity distribution of 5(A;, w) in the MI and SF phases as well as in the vicinity 
of the Kosterlitz-Thouless (KT) transition point, where the charge excitation gap closes. For 
large U, the spectral weight is mainly concentrated around u = U in the region 7r/2 < < vr, in 
agreement with the strong-coupling prediction. As U weakens in the MI phase, the distribution 
of the spectral weight broadens. At the same time, the maximum value of S{k, u) acquires a 
sizable dispersion (see upper-row panels). As the system reaches the MI-SF transition point we 
observe a significant redistribution of spectral weight to lower k values and, most notably, the 
excitation gap closes (see lower panels). In our previous work, we evaluated the scaling of the 
Tomonaga-Luttinger liquid parameter and determined the KT transition point to be located at 
t/U = 0.305(1) [3]. In the SF phase spectra for A; > (cf. the panel for t/U = 0.4), we observe 
an almost linear dispersion of the S{k,u}), which — in accordance with bosonization [11] and with 
Bogoliubov theory [12] — can be taken as a signature of the Bose condensation process. Our 
ID DDMRG BHM data are unsuggestive of two distinct (gapless sound and massive) modes in 
the SF phase. These two phases are found in mean-field theory and may appear for dimensions 
d> 2 where a true condensate exists in the SF ground state. 
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Figure 3. Intensity of the dynamical structure factor S{k, uj) in the k-uj plane for different t/U. 
DDMRG data were obtained for L = 32 system, using r] = 2t {ij = t) in the MI phase and at the 
KT transition point (in the SF phase). The red crosses mark the positions of the maximum in 
each k = 27rj'/L-sector. The green line for t/U = 0.4 marks the field-theory/Bogoliubov result 
uj{k) = y^eo{k)[eo{k) + 2U] with eo{k) = -2t{cosk - 1). 

To summarise, we have determined the dynamical structure factor S{k,uj) for the ID BHM 
with particle density p = N/L = 1 by means of unbiased numerical DDMRG calculations. 
As discussed for photoemission spectra previously [3], S{k,uj) agrees with the first-order 
perturbation theory result in the Mott insulator phase for U ^ t. Naturally, as the Coulomb 
interaction is lowered, noticeable deviations appear between both approaches, in particular the 
DDMRG S{k,u}) becomes dispersive and we find a substantial redistribution of spectral weight 
into the small A;-sector. In this regime, higher-order corrections have to be taken into account 
in our analytical treatment of the BHM. Approaching the SF state, the charge excitations gap 
closes and the maximum in S{k,u}) exhibits a linear dispersion. The quantum phase transition 
between MI and SF phases is located at about t/U ~ 0.3 and found to be of KT type. 
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